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Abstract 

Zonal flows are recognised to play a crucial role for magnetised plasma confinement. The genesis of 
these flows out of turbulent fluctuations is therefore of significant interest. We investigate the relative 
importance of zonal flow generation mechanisms via the Reynolds stress, Maxwell stress, and geodesic 
acoustic mode (GAM) transfer in drift-Alfven turbulence. By means of numerical computations we quantify 
the energy transfer into zonal flows owing to each of these effects. The importance of the three driving 
ingredients in electrostatic and electromagnetic turbulence for conditions relevant to the edge of fusion 
devices is revealed for a broad range of parameters. The Reynolds stress is found to provide a flow drive, 
while the electromagnetic Maxwell stress is in the cases considered a sink for the flow energy. In the limit 
of high plasma beta, where electromagnetic effects and Alfven dynamics are important, the Maxwell stress 
is found to cancel the Reynolds stress to a high degree. The geodesic oscillations, related to equilibrium 
pressure profile modifications due to poloidally asymmetric transport, can act as both sinks as drive terms, 
depending on the parameter regime. For high beta cases the GAMs are the main drive of the flow. This is 
also reflected in the frequency dependence of the flow, showing a distinct peak at the GAM frequency in 
that regime. 

PACS numbers: 52.25.Gj, 52.35.Ra, 52.65.Kj 
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I. INTRODUCTION 

Since the discovery of the H-mode 13] in magnetically confined plasmas a multitude of mech- 
anisms for the generation of the shear flow connected to the LH-transition have been proposed. 
They include amongst others ion-orbit loss effects, neoclassical effects, and turbulent flow gener- 
ation lO, 3I ^. Here we focus on turbulence as a source of shear flow generation. It was already 
early recognized that turbulence can lead to spontaneous self-organization of turbulent energy into 
sheared poloidal flows which in turn could reduce the transport significantly A conclusive 
computational demonstration of shear flow generation by turbulence in realistic geometry of fu- 
sion devices, which is sufficient for achievement of the LH-transition, has, however, not yet been 
achieved. In electrostatic turbulence the Reynolds stress is the main source of interaction between 
large scale flows and small scale turbulence. The Reynolds stress designates the radial flux of 
poloidal momentum, and a finite radial gradient of it will be an indication for a local condensation 
of momentum into a poloidal flow. In electromagnetic turbulence an additional source of poloidal 
flow generation has to be accounted for: the Maxwell stress, which arises from parallel momen- 
tum transport along perturbed magnetic field lines. Measurements of the Reynolds stress and its 
radial variation have been performed in several fusion devices with the purpose to identify it as a 
source of sheared poloidal rotation 1^. Recently, also the Maxwell stress, respectively, magnetic 
fluctuations and their cross-correlations have been measured in Reversed Field Pinch (RFP) |Q] 
and Tokamak [8] configurations. These measurements indicate that the Maxwell stress acts as a 
sink for poloidal flow energy. Finally, in the presence of toroidal magnetic field inhomogeneity 
the geodesic acoustic modes (GAMs) ^9^ interact with the poloidal flows in the system. In such 
cases the zonal flows show a residual oscillation at the GAM frequency. 

The purpose of this paper is to investigate these three different transfer mechanisms for zonal flow 
generation over a wide range of parameters as neither their strength nor their detailed (driving or 
damping) effect on the flows are a priori sufficiently clear. While the Reynolds stress is most often 
identified as a flow drive, there is considerable confusion about the role of the GAMs llfli [lli. 
The Maxwell stress is in low (3 situations rather weak, but it has been found to drain energy from 
the flow B and in high (3 situations it should ideally cancel the Reynolds stress [l^. 
This paper is organized as follows: In the following Section |ll| we present the turbulence model 
used for the computations. The next Section |llll is devoted to discussion of the various transfer 
mechanisms of energy between turbulence and flow motion in a low and a high beta case. We then 
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present global scalings of the transfer terms with coUisionality and plasma beta in Sec.|lVl Finally 
we discuss our results in the concluding section. 



II. ELECTROMAGNETIC TURBULENCE MODEL 

We investigate the detailed balance of drive and sink terms for global poloidal flows in a model 
for plasma turbulence in the edge region of magnetic confinement devices. Considering both elec- 
trostatic and electromagnetic effects, together with toroidal geometry and magnetic field curvature 
in a flux tube model, allows us to investigate the different turbulent momentum transfer terms re- 
sponsible for flow generation. 

The fluid equations for drift-Alfven turbulence in 3-dimensional flux tube geometry result from 
standard ordering based upon the slowness of the dynamics compared to the ion gyro frequency 
Q.i = eB/Mj and the smallness of the drift scale p.^ compared to the background pressure gradient 
scale length L±. These quantities and the sound speed Cs are defined by 

a = ^, c? = ^, p.. = |, L^ = \Vlogp,\-\ (1) 

where subscripts e, i refer to electrons or ions respectively, and the temperature is given in units 
of energy. Normalization is in terms of scaled dependent variables (electrostatic potential e<^/Te, 
electron density n/noo, parallel ion velocity u/cs, parallel electric current J/nooeCs). In addition 
the dependent quantities are scaled with the small drift parameter 5 = ps/L±, so that mainly terms 
of order one appear in the normalised set of equations. 

The scale perpendicular to the magnetic field is in units of p^; the parallel scale is L | = qR, with 
R the toroidal major radius and q the safety factor and the closed flux surface connection length 



27rL|| . The time scale is L^/cs- Further details for this system and geometry are given in Ref. lll4ll . 
In the present paper we use a constant background temperature for electrons, and assume cold 
ions. The quantity hqq is a normalizing density, while nQ{x) is the equilibrium plasma density 
having a finite gradient. In normalized units the radial profile of the density is dx\ognQ{x) = — 1. 
Thus X serves as the radial coordinate. Relative to the background magnetic field B the other 
perpendicular coordinate is y. Finally the parallel coordinate is denoted by z- 
As singly charged ions and quasi-neutral dynamics are assumed, uq and n refer to both the electron 
and ion density. It is important to note that hq and n equivalently describe electron density or 
pressure; we neglect temperature dynamics in this qualitative study because of the similarity in 
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physical character between the electron temperature and the "non-adiabatic" part of the electron 
density lisi . The model is described by the temporal evolution of the electrostatic potential ((|)), 
density perturbations (n), parallel current (J), and parallel ion velocity (u). Auxiliary variables are 
the vorticity (Q.) and the parallel component of the magnetic vector potential (A||): 

-^ + VE-Vn = !K{n)+V^\J+^ayl^, (2) 
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— + VE-y{no + n) = ^{n-(^)+V»{J-u)+^n^ln, (3) 



|- (^pA|| +jjJ^ +pvE • V7 = V|| (no + n - (^) - C7 , 



(4) 



e[^ + VE-Vu] =-Vn{no + n) , (5) 



dt 

with the vorticity and current (Ampere's law) J given by 

Q = Vl(^, 7=_v2A||. (6) 

The advective and parallel derivatives carry non-linearities entering through ^ and Ay, which - 
due to the description of the geometry - can be expressed in terms of a Poisson bracket 

{f g} = ——-—— 0) 

' dx dy dy dx 

in the ;cy-plane as 

Vi--V = a-}; V||=^-{pA||,-}. (8) 

The curvature operator Ki is for simple circular toroidal geometry written as 

Ki = -(Ob f sinz^ + cosz.^ ) , (9) 



dx dy 

and originates from compressibility terms of the form V • ( 1 /B^)B x V. Note that z takes values in 
the range [— TT : tt] and that the outboard mid-plane is located at z = 0. The perpendicular Laplacian 
is in the locally shifted metric \\\ written as 



and is due to —V ■ {b^B x 5 x , thus hiding magnetic shear in the shifting procedure. The 
viscous/diffusive terms ^ na-,l^n in Eqs. © and Q are introduced to provide sub-grid dissipation 
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of small scale dynamics. 

The parameters in the equations reflect the competition between parallel and perpendicular dy- 
namics, governed by the scale ratio £ = {qR/L^_)^. The electron parallel dynamics is controlled 
by 

|3=-^e, ^=-E, C = 0.51— ^ = v^, (11) 

where Xg is the electron collision time and the factor 0.51 reflects the parallel resistivity 1II6I1 . The 
competition between these three parameters, representing magnetic induction, electron inertia, and 
resistive relaxation, determines the response of J to the static force imbalance in Eq. Q. Due to 
the presence of in Eq. Q this adiabatic response has different character in different parts of 
the spectrum. The last physical parameter is (Og in Eq. © , reflecting the effects of magnetic curva- 
ture (equivalently magnetic gradient, in a toroidal model). An important note is that all magnetic 
induction dA^^/dt and flutter P{A||, ■} effects enter through the finite beta |3 = l^QPe/B^ or cl/v\, 
where va is the Alfven velocity, and P = (3e£. 

The density equation is augmented by damping layers in the left and right 5% of the radial do- 
main, regulating the poloidally averaged density, e.g. the profile modification, back to zero. This 
feedback control of the profile arranges for the average density profile to stay close to the one 
characterized by the originally defined gradient. 



III. ENERGETICS AND EVOLUTION OF FLOWS 



The equation determining the evolution of zonal flows is found from the vorticity equation 
Eq. (O by averaging over a flux surface as 

+ I (v.v,) - P ^ {B.By) + (0. {n sinz) = ^ , (12) 

where (•) = (l/lnLy) J%dz Jq^ dy- denotes the flux surface average. The E xB velocity is given 
by ve = (vx, Vy, 0) = (— dx^, 0) and the electric field connected to the poloidal flow is described 
by the potential 4>o(x) = ((])). 

Consequently Vo{x) = {vy) = dx^o and (Q.) = dxVo. The first contribution to the flow evolution 
is the Reynolds stress, which is the radial transport of poloidal momentum by radial velocity fluc- 
tuations. It demands a correlation between the two components of the fluctuating velocity, which 
may be due to a seed flow or a background gradient, as is the case here. The second contribu- 
tion arises from magnetic flutter. It can be interpreted as parallel current flowing radially along 



perturbed magnetic field lines. The third term is the acceleration of the flow due to interaction 
with density sidebands via the compressibility of the diamagnetic drift, associated with geodesic 
acoustic modes. Finally, viscosity on the right hand side of Eq. [T^l introduces a damping of the 
flow profile. 

We are interested in the energetics of the flow evolution. To find the evolution of the energy in the 
mean flow we multiply the vorticity equation Eq. Q by the flow velocity and integrate over the 
whole volume. We then obtain for the time evolution of the mean flow energy U := (1/2) / JxVq : 

dU 

— = 'K.+M +g +V , (13) 
with the quantities %_,M ,g , and V defined as follows: From the convection we find 

11= j Jx$o(vE-Via) = j dx{v^Vy)d_,Vo. (14) 

which is the Reynolds stress contribution to the flow drive. Correspondingly the parallel current 
and magnetic fluctuations lead to 

M =J dx^oi^wJ) = -P/ dx{B,By)d, Vo , (15) 

which is the Maxwell stress governing the energy exchange of the flow with magnetic fluctuations, 
where B^^ = 3yA|| and By = —dj^A^^. While the contribution of the normal curvature vanishes, the 
geodesic curvature results in a term associated with the geodesic acoustic modes (GAMs): 

g= J dx<i>o{(iiBsmsdxn) = —(Hb J dx{nVosms). (16) 

The coUisional damping finally is always a sink: 

V=-naj Jx4>o(v2 ^l) = j dx{d,VQf. (17) 

Flow generation by Reynolds stresses is well known to result from an average phase correlation 
between the velocity fluctuations in the drift plane spanned by the x and y coordinate axes. The 
tendency of convective structures to be tilted with a seed sheared flow makes the transfer term ^ 



generally positive, draining energy from the fluctuating motions to the zonal flows | |17|1 . 

It is worthwhile to note that in pure MHD turbulence there is an approximate balance between 
Maxwell and Reynolds stress [13]. From a local linear analysis of modes, neglecting the toroidicity 
of the equilibrium magnetic field, we obtain the following functional relationship between the 
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fluctuations in magnetic potential and electrostatic potential: 

{())Bky)/{h\k])+c 



, (18) 



" [a)fiV(P-/>^i)]/[^tl^i] + l 
with c = (i)/k\\. The dispersion relation has several branches (see Scott lisl '). In the limit of high 
P and neglecting effects of curvature, the Alfven branch of the dispersion relation dominates and 
c can be approximated by the Alfven speed va = P^^/^: 

Ai|=(^/VP- (19) 

As a consequence the Maxwell and Reynolds stress cancel in that regime, which is expressing the 
fact that Alfven waves do not transport poloidal momentum. 

From the plasma continuity equation ^ we find the evolution of the density sidebands, 

Tr-(nsmz) + Tr-(smzn— ) +a)5(sm ZTr-)=C0B(sm z^— )-(smz— ). (20) 
ot ox ay ax ax dz 

The contribution from the flow Vq = d(^o/dx in the first term on the right hand side of Eq. (l20t . 
describing the up-down asymmetric plasma compression due to poloidal rotation, couples with 
the zonal flow equation (fT^ and results in geodesic acoustic modes (GAMs) at frequency CO5 / \/2 
(Refs. IIHQQ])- 

Other terms in Eq. (l20b . along with coupling to the ion flow sidebands, 
may cause an acceleration of zonal flows in the presence of poloidally asymmetric particle fluxes. 
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known as Stringer- Winsor spin-up Ill0lll9l|20l|21|]. In this connection we also note that the energy 
transfer due to toroidal geometry into the energy of the fluctuating motions 

!K = ydxv\ (21) 

is given by 

— J dx^Kiln) = —(Ob J dx I sin^n — -|-cos5n — 

This indeed indicates the tendency towards a ballooning structure of the fluctuations, since this 
term drives velocity fluctuations when the turbulent plasma transport is radially outwards from 
the torus axis and poloidally towards the out-board mid-plane. This geodesic transfer process was 
recently revisited in Ref. [11], where it was claimed that the GAM transfer is generally from the 
zonal flows through the density side-bands to the turbulent fluctuations. 
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IV. STRUCTURE OF ELECTROMAGNETIC FLUCTUATIONS 



To address the simultaneous action of the energetic transfer effects we resort to three- 
dimensional numerical computations of the four- field model eqs. ©- (|51) on a grid of usually 
64 X 256 X 32 points, with dimensions 64 x 256 x 2n in x, y and z, respectively. Some runs were 
repeated at higher resolution 128 x 512 x 32 to ensure convergence. The numerical scheme uses 
a symm etry , energy and vorticity conserving discretisation of the bracket structure of the nonlin- 
earities 11221] with the curvature terms cast into bracket form as well. Time stepping is performed 
using an explicit third order stiffly-stable scheme |23|], with viscous terms treated implicitly using 
operator splitting. For more details on the numerical implementation see [24]. 

Nominal parameter values typical for tokamak edge plasmas are £ = 18750, /i = 5, s = 1, 
(Og = 0.05 and iuq, = = 0.025. 

For parameter scans we varied P in a range between 0.1 - 30 and v from 0.5 - 7.5. The 
scaling with |3 is especially interesting, as the drift- Alfven system has the property that the nature 
of the turbulence changes with the degree to which the system is electromagnetic. This feature 
was demonstrated numerically by Scott lll8l] and Naulin ll24ll and also experimentally by Lechte 
et al [25]: The transition manifests itself in a change of the phase relationship between density 
and potential fluctuations, which varies for low values of ky from a small phase angle in pure drift 
wave dynamics to 7r/2 in the MHD drift-ballooning regime. This is exemplified in Figure[Tl which 
shows the phase probability distribution function as function of poloidal wavenumber for the cases 
P = 0. 1, V = 2.295 and P = 30, v = 0.5. While in the low P case, the phase angle is always small, 
for the large P cases we observe a much broader phase relationship and a generally larger phase 
angle. The regime of dominating MHD ballooning instability is first reached at for the edge very 
high values of P > 30 

A time-trace of the kinetic energy, ac, of the fluctuating motions and the zonal flow energy 
U = J dx{l/2)VQ is presented in Fig.|2l It is seen that while saturates after about 100 time 
units, the saturation of the flow takes place much later. Thus, all computations were run to times 
t = 5000, with time averages taken in the interval from time t = 1000 to the end of the simulation, 
to ensure a statistical steady state of fluctuating quantities. Moreover, from Figure|2lit is observed 
that the energy in the zonal flows is only a fraction of the total kinetic energy. This underlines 
the fact that no strong self-organized H-mode like transport barriers are formed in this system. 
However, zonal flows do form and influence the profile of transport as well as the density profile. 
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Due to the change in turbulence character connected to P we will now proceed and present in more 
detail two runs, the low beta P = 0. 1, v = 2.295 and a high beta P = 30, V = 0.5 case. 

Figure |3l shows a gray-scale plot of the zonal flow profile Vo(x, t) and the zonal density (n) (x, t) 
in time, where in both cases we omit the damping layers in the plot. It is clearly seen that the zonal 
flows are radially localized and while exhibiting some fluctuation features, the flow profile is rather 
persistent in time. The zonal density shows some imprint of the zonal flow in terms of slightly 
elevated density levels in the vicinity of high flow shear, but fluctuations in the zonal density are 
more pronounced. It is worthwhile to remark that at about t ^ 3750 a weakening of zonal flows is 
clearly noticeable and prompts a transient radially propagating feature in the zonal density. This 
provides us a visualisation of the interplay between flows, transport and the density profile. 

In Fig. m we present time traces of the energy transfer terms into the zonal flow and the rate of 
change of the zonal flow, together with the numerical error obtained by comparing the sum of the 
energy transfer terms with the computed actual change rate of the flow: 

5F(0 = ^ +g + V) . (22) 

dt Num. 



Here dVo/dtNum. is evaluated to second order in time i26\ . The error 8F{t) is of the same order, 
and for P = 0. 1 the Maxwell stress energy transfer term fAf , which is negligible for this low 
value of P, when compared to each other energy transfer term, is of the same size. The statistical 
nature of the fluctuating flow drive terms is observed, as the balance between the transfer terms 
is only reached on a long time average, whereas on an instantaneous view the transfer terms can 
deviate significantly from their means. Here viscous damping 'U and GAMs g serve as sinks 
for the flow energy, which is solely driven by the Reynolds stress. While both damping terms 
vary on a rather slow time scale, the Reynolds stress and with it the resulting rate of the flow 
change vary on the faster time scale of the turbulence. Figure |5] shows correspondingly time 
traces of selected fluctuating quantities obtained at a single point and of flux surface averaged 
quantities at the same radial position. While the fluctuating quantities are all varying on the fast 
scale, the flux surface averaged ones vary significantly slower. Zonal density and zonal magnetic 
potential show, however, variations on a time scale of about 250 and 80 time units, respectively. 
The zonal flow Vq shows some fast scale jitter, but varies only slowly over the shown 1500 time 
units. Thus, to investigate that time behavior in more detail we present in Fig. |6lthe frequency 
spectra of flux surface averaged quantities connected to GAMs, as the density and the parallel ion 
velocity together with quantities related to zonal flow dynamics. The zonal density (n) and the 



flux-surface averaged parallel ion velocity (m) show both a pronounced peak at a low frequency of 
about (0 ^ 0.025. This peak is clearly associated to the ideal geodesic oscillation around QOgam ~ 
(£)b/V2 = 0.035. It is worthwhile to note that the GAM frequency arises by combining Eq. (fT2|) 
with Eq. (1^ . and the ideal GAM frequency arises from the relation 

(UBisinh^) = ^(Ob{[1 -cos(2z)] ^) ^ ^(ObVq (23) 

if the flux surface average of the term {cos{2z)dx^) disappears exactly, as it would be expected 
for a fluctuations Vy being homogeneous along the parallel coordinate z. In toroidal geometry 
the <^{z) and thus Vy{z) = however, show in general a distinct ballooning feature, resulting 
in higher amplitudes around position z = than for z = ±71. We thus expect the GAM peaks in 
the spectra to be shifted from the ideal (Ugam, with the direction of the shift depending on the 
preferential direction of local flows Vy, and the width of the frequency shift depending on the 
ballooning properties of the velocity fluctuations Vy{z). For our present parameters and a probe 
location one third into the radial x domain we experience a downshift by an additional factor of 
approximately a/1/2. 

For the zonal flow Vq and the zonal vorticity, we observe that the zero frequency mode domi- 
nates the poloidal flow spectrum. At the frequency of the zonal density feature we observe even 
a small dip in the flow frequency spectrum, this supporting the observation that GAMs are a sink 
for the flow energy in that parameter regime. 

For the high |3 = 30 case the situation looks differently as is clearly seen from Fig.|7| The zonal 
flow profile is now broader and reveals much less persistence than in the low beta case (compare 
Fig. 121). Correspondingly the zonal density shows also a less pronounced radial structure and the 
characteristic time of the fluctuations seems to be of similar size for both zonal flow and zonal 
density profile. 

The energy transfer rates M , and g , shown in Fig. [S] reflect this change in behavior. We 
first observe that the Maxwell stress is now of finite size and a significant sink for the flow energy. 
It is very well correlated to the Reynolds stress in time, which still acts as a flow drive. For an ideal 
high beta MHD case in linear geometry without magnetic field curvature the balance between M 
and ^ is known to be exact with no preferred flow direction in the dynamics, as seen in Eq. (fT9b . 
Here the resulting energy transfer from Maxwell and Reynolds stress is close to zero and the 
resulting change in the effective zonal flow drive is dominated by the GAMs g , which now acts as 
a driving term. Consequently the resulting rate of change for the flow evolves mainly along with 
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the change in the GAM drive and shows only a minor additional variation on the timescale of the 
Reynolds and Maxwell stresses. 

The time traces of fluctuating and averaged quantities shown in Fig. 13 reflect that behavior. We 
first note that the fluctuating quantities now show a somewhat more pronounced slower frequency, 
which reflects the frequency observed in the time evolution of the zonal quantities. All fluctuations 
are larger by about a factor three compared to the low beta P = 0. 1 case. The zonal density is up 
by a factor two to threee and a slow oscillation is clearly observed in both the zonal density as 
the zonal magnetic potential. This slow frequency now is also found in the zonal flow time trace. 
These features get more obvious in the frequency spectra depicted in Fig. [lOl A pronounced 
low frequency behavior is now seen also in the flow related quantities, namely in Vq, which here 
exceeds the zero frequency component by about a factor two. In these situations the flow is not 
stationary (zero frequency) but is, compared to the turbulence, a slowly varying structure. The 
slow frequency of the flow is close to the ideal GAMs frequency QOgam, revealing the flow drive 
by this process. These results are generally in agreement with experimental observations, that 
show a modulation of the zonal flows at the frequency of the GAM oscillation 



V. PLASMA BETA AND COLLISIONALITY SCALINGS 



Here we present results concerned with the scaling of the different transfer terms with colli- 
sionality and plasma beta. In Fig. [TT] we show the three main transfer terms as a function of |3 
for a low coUisionality of v = 0.5. For increasing P the Reynolds stress drive gets slightly 
weaker, but is in all cases a drive. This indicates that the described flows do not decay through a 
Kelvin-Helmholtz like instability mechanism, which would make the Reynolds stress a sink term. 
The Maxwell stress 5Vf starts close to zero and is always a sink term. It grows as expected in sig- 
nificance with increasing |3, and for large P is the dominating sink for the flow energy. The GAM 
transfer, g , starts out as a sink for the flow energy at low beta, but with increasing beta it looses 
its importance as a sink. Finally, g becomes positive for the high beta P = 30 case, e. g. the GAM 
acts as a flow drive. 

We then look at the scaling of flows and energy transfers with coUisionality in the two cases 
of low and high beta. The results are presented in Fig. [121 We find that the saturation level of 
zonal flow energy decreases with rising v for low P and is in general by an order of magnitude 
smaller than the fluctuating kinetic energy K = j '^x(l/2)(V^^)^. The system is mainly gov- 
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emed by vortex dynamics 112811 . where the zonal flow contribution plays an important part for 
self -regulation of the fluctuation amplitudes, but is not dominating the energetics to an extent that 
it would completely suppress the turbulence. Further, we observe that the fluctuation energies ^ 
and ^ = J dx{l/2)fP' and the turbulent particle transport r„ = / dw^n both increase with paral- 
lel resistivity v. This stronger turbulence level is due to the increasingly non-adiabatic nature of 
the electron response due to collisions, which increases the nonlinear drive of drift modes. The 
energy transfer terms behave as follows: The Reynolds stress decreases with increasing coUision- 
ality, accordingly with the decaying flow energy. The Maxwell stress energy transfer M is always 
negligible at low beta and thus the flow energy is dissipated through the viscous terms and the 
geodesic channel. 

Energies , and the transport r„ approximately double when changing to |3 = 30 and taking 
the system from the drift into the ballooning regime. The magnetic flutter effect is then an impor- 
tant cause for the non-adiabatic response on the electrons. The influence of v is accordingly much 
smaller and results in less variation for this high (3 case. 

We observe that the Reynolds stress is always a drive (5?. > 0) for the flow. At high beta 
the Maxwell stress is important and the balance between Maxwell and Reynolds stress becomes 
obvious as 1 51/ I ^ is reached and exceeded. Finally the electromagnetic flutter effect is taking 
out more energy from flows than is injected by Reynolds stress spin-up. The flow sustained in 
this high beta regime is now maintained by geodesic transfer into the flow. The driving effect 
by GAM oscillations on the flow is more pronounced for higher v and higher levels of transport: 
GAM transfer is closely linked to the energy of density fluctuations in the m = ±1 sideband 
(and thus to energy in all other scales that couple by three-wave interaction to this sideband), 
which directly scales the transfer term g . For higher resistivity the relative importance of g is 
thus enhanced in the same amount as n increases both due to a more resistive as well as more 
electromagnetic electron response. The strong magnetic flutter B± in high beta turbulence can 
cause significant chaotic deviations of field lines radially and poloidally from the flux coordinates 
that were defined for an unperturbed magnetic field. The parallel coupling is thus able to connect 
regions of neighboring radial domains where the amount of radial overlap is rising with (3. Zonal 
structures may be efficiently destroyed for a strong flutter effect, and the aligning character of 
Reynolds stress on zonal flows is counteracted, thus also the radial structure of the zonal flows is 
less pronounced in the high beta case, see Fig. |7l The drive of vortices on the drift wave scales is 
of course still maintained by the free energy in the background density gradient, and the energy 
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on drift scale density fluctuations is even increased by the destabilising magnetic flutter effect 
on the non-adiabatic parallel dynamics. The cascade in density structures is generally a three- 
wave interaction that is on statistical grounds essentially down to smaller scales, but is by more 
infrequent events also able to feed scales in the m = ±1 geodesic sidebands and the m = zonal 
mode. The geodesic transfer pathway is thus open in both directions: a strong drive of the GAM 
mode by zonal flows for low beta on the average drains energy out of the flow to smaller scales. 
For high beta the GAM energy is, however, converted to a certain extend into ^ oscillations that 
supply the zonal flows. 

VI. CONCLUSION 

We have performed a detailed investigation of the zonal flow drive in drift- Alfvn turbulence for 
parameters relevant to the edge region of hot plasmas in toroidal devices. We have identified three 
main mechanisms for the interaction of the zonal flows with the turbulent fluctuations; namely 
the electrostatic Reynolds stress, the electromagnetic Maxwell stress, and the geodesic acoustic 
mode, GAM, coupling. The main results are summarized as follows: For low beta cases the 
Maxwell stress is negligible and the Reynolds stress is the only driving term of the flow, whereas 
the GAM coupling provides a sink for the flow in addition to the viscousity. For the case of high- 
beta plasmas, however, the Maxwell stress becomes significant. It acts as a sink for all the cases 
we have investigated, and it efficiently cancels the driving effect of the Reynolds stress. In this 
parameter regime the flow is mainly sustained by the GAM coupUng, that now acts as a drive 
opposing the viscous damping. 

We should therefore emphasize that from an experimental point of view, measuring Reynolds 
stress exclusively as an indication for flow generation is in general not sufficient. The electro- 
magnetic Maxwell stress is important already at a moderate edge beta parameter, and will be even 
more important for ITER like plasmas with higher edge [3 at reduced coUisionality. This clearly 
opens a demand for additional measurements of the Maxwell stress. 

There is a clear trend in the computational results that assign more importance to the GAM 
oscillation at a high level of transport. The GAMs present a driving mechanism for the flows if the 
transport is sufficiently inhomogeneous: with an increasing ballooning character of the turbulence 
the GAMs are further excited and can ultimately drive flows. 

The frequency spectra of the zonal flow clearly show a dip or a peak in the (Ugam frequency 



13 



range, depending on the sink or drive role of the GAMs for the flow evolution. Measurements 
of the frequency spectrum of the zonal flow should thus be able to distinguish between these two 
scenarios and provide further insight into the importance of GAMs for the flow, and finally for 
H-mode formation. 

Finally we note that our numerical results are for the high beta case partially in disagreement 



with recent results by B. Scott llli regarding the specific role of the Maxwell stress Il29l] 
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FIG. 1 : Phase angle probability between density and potential fluctuations for the low P = 0. 1 (top) and 
high P = 30 (bottom) case. 
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FIG. 2: Kinetic energy and energy in zonal-flow component U over time for the low P = 0.1 case. 
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FIG. 4: Energy exchange terms, flow change rate, and numerical error for (3 = 0.1. Reynolds stress is the 
key drive and GAMs are acting as a sink for the flow energy. 
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FIG. 5: Fluctuating quantities (top) and fluxsurface averaged quantities (bottom) for P = 0.1, measured at 
X = Lx/3 and on the outboard midplane. 



21 



0,3 



0,2 



0,1- 



I I I I I Mil 1 I I I I 



I I I I I 



p = o.i 



(n)(L,/3) 
1000(w)(L^/3; 



0,0' ' 

10 10 



0,08 




10"^ 10"^ 
coL /c 

n s 



10 



0,06- 



0,04 



0,02 



0,00 



I I 1 1 1 



I I 1 1 1 



I II 1 1 1 



I I 1 1 1 



TT 



{n){L,/3) 

0.21^(L^/3) 




coL /c 

n s 



FIG. 6: Frequency spectra of quantities associated with GAM oscillation (top) and flows (bottom) for 
P = 0.1. The vertical Une indicates the ideal GAM frequency cogam- 
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FIG. 8: Energy exchange terms, flow change rate, and numerical error for (3 = 30, showing the dominating 
influence of the GAMs for flow drive in that regime. 
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FIG. 9: Fluctuating quantities and fluxsurface averaged quantities for P = 30, measured atx = Lx/3 and on 
the outboard midplane. 
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FIG. 10: Frequency spectra of quantities associated with GAM oscillation (top) and flows (bottom) for 
P = 30. The vertical line indicates the ideal GAM frequency (Ogam- 
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FIG. 12: Scan over collisionality v for low p = 1.0 (top) and high P = 30 (bottom). The left side shows 
energy in the density fluctuations P, kinetic energy K, flow energy U and particle flux F,,. The right side 
depicts energy transfer terms. 
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